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Abstract 

For autonomous systems it is well known how to extract tunneling probabilities from wavepacket 
calculations. Here we present a corresponding approach for periodically time-dependent Hamilto- 
nians, valid at all frequencies, field strengths, and transition orders. After mapping the periodically 
driven system onto a time-independent one with an additional degree of freedom, use is made of the 
correlation function formulation of scattering [J. Chem. Phys. 98, 3884 (1993)]. The formalism is 
then applied to study the transmission properties of a resonant tunneling double barrier structure 
under the influence of a sinusoidal laser field, revealing an unexpected antiresonance in the zero 
photon transition for large field strengths. 

PACS numbers: 03.65.Xp, 05.60.Gg, 72.40.+W, 73.40.Gk 
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Recently, the quest for experimental controllability of transport processes on the 
nanoscale has triggered a number of theoretical studies on the influence of periodic external 
fields on molecular scale conductors \u,M as well as on semiconductor resonant tunneling het- 
e..os..uc.u.es SQQ a„d ,ua,.u,„ dot. Q. Diffe.e.t theo.e.ca. tecKni^e. .a^g.ng .o-n 
Green's function approaches 

propagation methods have been employed to calculate either transmission probabilities 
or directly the currents for the considered setups. 

In the seminal work of Biittiker and Landauer j3| an analytical approach to calculate the 
transmission probability of a sinusoidally driven potential barrier was made. The perturba- 
tive analysis allowed the calculation of central and side band transmission probabilities at 
the incoming energy and at the incoming energy plus or minus the photon energy. An anal- 
ogous perturbative way of extracting those "driven" transmission probabilities was shortly 
afterwards used in numerical ^ and analytical P| studies of resonant tunneling. The non- 
perturbative calculation of tunneling probabilities under absorption or emission of photons 
is a necessity if the influence of strong fields on the tunneling process is to be investigated, 
however. These probabilities play a central role in extending the Landauer formalism to 
conductance calculations in the presence of driving . A viable generalization of standard 
approaches to make the determination of tunneling probabilities for arbitrary field strengths 
possible is therefore needed. 

A very intuitive way to calculate the scattering S{E) matrix as a function of energy, from 
which the tunneling probability T{E) can be extracted, has been given in the framework 
of reactive scattering theory for autonomous systems This approach relies on the 

solution of the time-dependent Schrodinger equation and the subsequent Fourier transform of 
cross-correlation functions. According to the reactive scattering situation one distinguishes 
between internal and translational degrees of freedom, where the translational one usually 
represents a tunneling degree of freedom. In the following, this approach will now be modified 
to include, instead of the internal coordinate a new coordinate representing time in the case of 
a periodically time-dependent system. The absorption or emission of a quantum of radiation 
can then be accounted for by calculating the appropriate cross-correlation. We will elucidate 
the formalism by applying it to a resonant tunneling Hamiltonian modeling e. g. a weak link 
metal- metal contact. 

In order to investigate the phenomenon of tunneling under the influence of an external 
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periodic field we follow an S-matrix approach. Therefore, the Hamiltonian of the driven 
system with a periodic driving term of the form f{x) sm{ujt) has to be made explicitly 
independent of time. This can e. g. be achieved by defining a new pair of canonical variables, 
Q = cot and pe, with the frequency cu of the external field jl2|. The autonomous Hamiltonian 
for a driven particle of mass m then depends also on the new momentum and coordinate 
and reads 
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H{p,, x,pe, 0) = ^ + Vi^) + fi^) sin(e) + upe- (1) 

This Hamiltonian can easily be quantized by the usual replacements p^^e = —iTT^dx,®- In the 
following, we assume that this has been done by keeping in mind the operator nature of the 
Hamiltonian (in a different notation also known as the Hamiltonian 

If the coupling term and the potential V{x) are localized, i. e. V"(x), f\x) ^ for |x| — * cxo, 
then the Hamiltonian at large values of |x| can be written as a sum of a kinetic energy term 
plus a term independent of x 
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H ^ i/ofe, x,P0, 0) = ^ + Kpb), (2) 

where the eigenf unctions of h{pQ) = ujpQ are exponentials, in their normalized form given 
by 

Xn(e) = ^e^"®, (3) 

and living on a finite support < < 2tt. In the following we will refer to n as the number 
of quanta in the external field, due to the energy eigenvalue En = nhu corresponding to the 
eigenfunction above. The eigenfunctions of the translational part of the Hamiltonian are 
energy normalized plane wave states with wavevector k. For the asymptotic Hamiltonian 
Hq this yields the direct product eigenstate and energy 

^n,E = Xn(e) J — -=e*^^ E = nhu + h^e/{2m). (4) 

The eigenstates of the full Hamiltonian H are created by applying the M0ller operator, 
defined by fl± = limj^::poo exp{ii7t//i} exp{— ii^oV^}) to the asymptotic eigenstate 

^n,i? = ^±'^n,E (5) 

and serve to define the on-shell S-matrix by 

SK,n';L,niE)S{E ' E') = (6) 
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with a slight generahzation of notation by including the "channel" indices L,R, denoting 
to which side of the barrier the eigenstates correspond. As noted in eigenstates can 
be represented by Fourier transform of propagated wavepackets. This idea applied to Eq. 
(jni) has been used later-on to extract the scattering matrix using wavepackets ^|. This 
reasoning can now be taken over to the present situation with the only difference that the 
eigenstates corresponding to the internal motion are states of the field variable G defined 
in Eq. Q. The wavefunction to be propagated is a Gaussian wavepacket centered in the 
asymptotic regime of the translational degree of freedom times an eigenfunction in the field 
variable 

0L,n(a;, O) = gh{x)xn{Q)= ( — ) 

exp{ - a{x - Xaf + ipa{x - Xa)}Xn{Q)- (7) 

For the S-matrix this leads to the expression 
where a cross correlation function of the form 

CR,n';L,n(^) = (0R,n' I exp{ -zi^t//l} 1 0L,n) , (9) 

is used, with a final state wavepacket defined analogously to Eq. ((Tj) but to the right of the 
barrier. Furthermore, the arbitrariness of the initial and final wavefunctions is removed from 
..eS-n,atH.exp.e.io„ by t.e„o.n,a,.a«o„ facto. 



VkAE) = / e-^^-^g^{x)dx, (10) 



with kn = \j2m{E — En)/h and a corresponding formula for rji^y. We stress that one does 
not need to use M0ller states as initial and final wavefunctions as long as they are located 
far in the asymptotic regions of the x-coordinate, such that the M0ller operator becomes 
the unit operator. Furthermore, we note in passing, that an approach to extract scattering 
information for periodically driven quantum systems, similar in spirit, but based on Floquet 
theory has e. g. been pursued in jl^, 

In the following we calculate transmission probabilities of a periodically driven sys- 
tem by numerically solving a 2 degree of freedom time-dependent Schrodinger equation 
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with Hamiltonian (^J). The x degree of freedom shall correspond to a tunneling elec- 
tron coupled to the external field, represented by the B degree of freedom. Without 
driving, the electron is supposed to move in the symmetric resonant tunneling potential 
V{x) = [Vn(l + exT)^ —l3xr})] /\l + exp|/3(|x| — Xc)}] — Kes exp{— 7x^} introduced by Bringer 
et al. |l7[ and applied in 171, |l8| to study inelastic tunneling in the presence of coupling to a 



harmonic oscillator mode. We are using the parameter values Vq = 10 eV, (3 = 4 a.u., 7 = 1 
a.u., Xc = 4: a.u. and Kes = 14 eV taken from j3] in the following. The undriven potential 
is then of the double barrier type, supporting a resonant tunneling "level" with a width of 
about 0.4eV at an energy oi E ^ 4.9eV, well below the barrier top. 

The periodic field is applied by using the function /(x) = xV^6(xc — localized in x 
around the double barrier potential. This model is motivated by a lead-molecule-lead setup 
irradiated by a laser field, the field-matter interaction being treated in dipole approximation 
and cut off due to the metallic nature of the leads. A dipole driving without cutoff could 
in principle be dealt with by applying a Kramers-Henneberger (KH) transformation to an 
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oscillatory reference frame 1^. This leads, however, to a complicated dependence on 
position and time of the potential V{x, t) and has to be dealt with by a transformation to the 
momentum gauge . The field parameters are chosen to make the expected effects clearly 
visible. For the frequency we thus choose hu = 1 eV, while the strength of the coupling 
to the external field will be varied in order to study its influence on the transmission 
probabilities. We will use parameters for which a perturbative approach becomes more and 
more questionable. 

Numerical results will be presented for the dynamics of wavefunction (0) with Gaussian 
width parameter a = 1 atomic units (a.u.), an initial center of position x^ = 10 a.u. and 
an initial center of momentum pa = 0.7 a.u., having a sufficiently large overlap with the 
transmission resonances. We start out by assuming that the G degree of freedom is initially 
in the n = state. Although, the number of quanta in the field has to be very large in 
order to justify the classical treatment of the driving term in Eq. (^, the total quantum 
wavefunction can be multiplied with a factor exp{i{m — n) [0 — ut] } and is still a solution of 
the time-dependent Schrodinger equation under the Hamiltonian (0). The exponential factor 
changes the number of quanta from n to m, justifying our "shift of the origin" to n = 0. The 
final states, which the time-evolved wavepacket is to be overlapped with, are again direct 
products as in Eq. (jZj), but located on the right side of the barrier (with x/3 = —Xa and 
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Pl3 = Pa), with eigenstates in the 9 degree of freedom labelled by n'. We stress that there is 
no restriction on the transition (or "side band") order n' imposed by our formalism. In the 
following, we focus on the cases n' = 0, 1,2, however. The grid employed in the numerics 
has 8192 points in the x direction for the electron and 32 points in the B direction. In order 
to avoid unphysical reflections from the grid boundaries we have employed an absorbing 
negative imaginary potential smoothly turned on at large values of 

In Fig. the correlation functions entering the Fourier transforms in Eq. (jH)) are depicted 
for increasing coupling strengths of Vc = 0.5, 1, 2 eV as a function of time in atomic units. 
In panel (a) for the 0-0 transition, it can be seen that the coupling to the field leads to a 
change in the long time tail of the oscillation of the correlation function cr_o:L,o('^) setting in 
at about 50 a.u.. For the cases n' = 1, 2, i. e. excitation of one, respectively two photons of 
the field during the scattering process, the oscillations in the time signals last longer than 
in the n' = case and die out at times of around 800 a.u. (not shown). The amplitude of 
the oscillations is substantially smaller than in panel (a). The overall shape of the signal in 
Fig. ^a) is not changed very much by the coupling. This does not imply that the resulting 
transmission will be unaffected, however. In a semiclassical study of undriven Id-tunneling 
it has e. g. been shown that, due to the effect of the normalization, minute differences 
between the semiclassical and the full quantum signal can lead to pronounced differences in 
the corresponding tunneling probabilities [l^. A similar effect of real, physical origin can 
be observed in the following due to the coupling to the light field. 

To this end, we now calculate the tunneling probabilities from the correlation functions 
of Fig. HI according to 

Tn'o{E) = \S^^^n';L,n=oiE)\'^. (11) 

as a function of (total) energy by using Eq. (jHl) for the S-matrix. Due to the coupling between 
the field and the electron, the final photon number n' may increase (decrease), accompanied 
by an equivalent energy loss (gain) of the electronic system, which has thus "emitted" 
( "absorbed" ) a corresponding number of photons. The numerical results are depicted in Fig. 
121 For the 0-0 transition, and for Vc = 0.5 eV, an isolated resonance with barely below unit 
transmission (solid line in panel(a)) at around 4.9 eV is observed, close to the unperturbed 
resonance energy. For increasing coupling strength, the peak in the transmission curve 
acquires a shoulder at about hu below the original peak and a corresponding dip to the right 
(panels (b) and (c)). The transmission probabilities with n' = 1,2 increase with increasing 
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field strength, exhibiting a doubly peaked structure, with the dominant peak shifted from 
the unperturbed resonance position by huj to higher energies and a shoulder blue-shifted by 
around 2huj which is most clearly visible for = 1 eV. The peak shifted by huj describes 
the scattering of an electron coming in with an energy higher than the resonance but by 
loosing one quantum of energy closely matching the resonance energy leading to an increased 
transmission probability. For = 2 eV, however, the resonances are already considerably 
broadened having almost lost the separated peak structure. Furthermore, the maximum of 
the 0-0 transition is shifted considerably to the left and the antiresonance character of the 
dip at around 6 eV is clearly emerging. This behaviour can no longer be described and 
understood perturbatively and is in sharp contrast to the case of inelastic coupling to a 
harmonic oscillator where there are no dips but only peaks at integer multiples of the 
oscillator frequency to the right of the unperturbed resonance. 

Finally, to check consistency, we have studied the case of a single field quantum being 
present initially by using n = 1. As a function of relative translational energy Et = E — nhuj 
in the incoming channel, the corresponding transmissions (not shown) then match exactly 
the ones shown in Fig. |21 i. e. Tu^Et) — * Tqq{E) and T2i{Et) — Tio{E). As has been 
noted before, the field eigenstate quantum number does not influence the scattering process. 
To compare the results directly, one has to consider the relative translational energy as the 
independent variable, however. Furthermore, in the present case, there is also an equivalence 
between e. g. Toi{E) and Tio(-E) (both as a function of absolute energy). The special form of 
the potential together with the driving term make the Hamiltonian obey generalized parity 
symmetry (see e. g. 20] for a recent discussion of this symmetry in a similar context). Thus 
it can be shown by using Eqs. (PfTTj) that the two transmissions have to be equal. We have 
checked that this is indeed the case also in our numerical results (Tqi not shown). 

We have shown that nonperturbative, exact numerical calculations of periodically driven 
tunneling probabilities, to any transition order, can be performed by employing an au- 
tonomous Hamiltonian in extended phase space. Within the wavepacket formalism of scat- 
tering pioneered in the probabilities for tunneling under absorption or emission of field 
quanta can thus be extracted from suitably chosen cross-correlation functions. The presented 
approach can be applied to all driven scattering problems, as long as the field coupling f{x) 
is localized in space. A way to deal with pure dipole driving in a scattering situation has 
been devised in jl^ by using a KH transformation which could also be employed in the 
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present formalism. 

For resonant tunneling, in sharp contrast to the case of inelastic coupling to a harmonic 
oscillator, the field case leads to an antiresonance in the 0-0 transition at high frequencies. 
Furthermore, in addition to the resonant tunneling problem, also systems like driven quan- 
tum dots j6| or molecular wire ratchets [Ij , and nonresonant barrier transmission problems, 
e. g. through semiconductor heterostructures, could be studied using the new methodology. 

The author is indebted to Sigmund Kohler for valuable discussions and would like to thank 
the "Deutsche Forschungsgemeinschaft" for financial support through the "Forschergruppe 
FOR 335", project C4. 
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FIG. 1: Real parts of correlation functions CR^n';L,n 

for Vc = 0.5 eV (full line), Vc = leV (dashed 
line), Vc = 2 eV (dotted line) as a function of time in a.u.: (a) n = n' = (b) n = 0,n' = 1 (c) 
n = 0, n' = 2 
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FIG. 2: Transmission probabilities T„/„(£') for Vc = 0.5 eV (a), Vc = I eV (b), Vc = 2 eV (c) as 
a function of energy (translational equals absolute in the present case) of the electronic degree of 
freedom in eV: n = n' = (full line), n = 0, n' = 1 (dashed line) n = 0, n' = 2 (dotted line). 
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